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Modified gravities are tlie most favorable candidates for explaining the current accelerating ex- 
pansion of the Universe. In this regard, we study the viability of an alternative gravitational theory, 
namely f{R, G), by imposing energy conditions. We consider two forms of f{R, G), commonly dis- 
cussed in the literature, which account for the stability of cosmological solutions. We construct the 
inequalities obtained by energy conditions and specifically apply the weak energy condition using 
the recent estimated values of the Hubble, deceleration, jerk and snap parameters to probe the 
viability of the above-mentioned forms of f{R, G). 
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I. INTRODUCTION 

As is well known in general relativity, the energy condi- 
tions are often required in the proofs of various important 
theorems about black holes, such as no hair theorem or 
laws of black hole thermodynamics [ij . The energy con- 
ditions are systematically obtained when one refers to the 
Raychaudhuri equation, in which the attractive charac- 
ter of gravity is reflected through the positivity condition 
R^uk^k'^ > 0, where i?^^ is the Ricci tensor and fc'' is 
any null vector. In general relativity, through the Ein- 
stein field equation, this condition on geometry is trans- 
lated to a corresponding condition on matter, namely null 
energy condition T^^k'^k'^ > 0, where T^i^ is the stress 
energy tensor. In particular, the weak energy condition 
(WEC) assumes that the local energy density is positive 
and states that T^i.C/'^C/'^ > 0, for all timelike vectors 
If^. For instance, in the case of a perfect fluid 



(1) 



we have p > and p+p > 0. By continuity, the WEC im- 
plies the nuU energy condition (NEC), Tf^^k^k"" > [J. 
The energy conditions have been widely studied in in the 
context of modifiedgravity, such as f{R), f{G), f{R,T) 
and /(T) gravity [J II]. The motivation for the study 
of modified theories of gravity is clear and strong. As 
we know, Einstein's theory of General Relativity (GR) 
can not explain the late-time accelerated expansion of 
the Universe, unless an unknown dark energy element 
is introduced in the field equations In the context 
of GR, the attempts to explain the observed speed-up 
of the Universe, have introduced some modifications of 
Einstein-Hilbert action by adopting a general function 
of the scalar curvature in the gravitational Lagrangian 
density as f{R) The motivation for this procedure 
consists of the analysis of strong gravitational fields near 
the curvature singularities and considering the consistent 
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candidates of a fundamental theory of quantum gravity. 
Indeed, string/M-theory predicts that scalar field cou- 
plings with the Gauss-Bonnet invariant G are important 
in the appearance of non-singular early time cosmologies. 
These motivations may also be considered in the context 
of late-time acceleration of the Universe 0-01 • 

Recently, a new generalized modified Gauss-Bonnet 
gravity, whose action contains a general function of R 
and G as f{R,G), has attracted considerable attention 
Q . Besides its stability, this attention is due to its ability 
to describe the present acceleration of the universe as 
well as the phantom divide line crossing and transition 
from acceleration to deceleration phases. Two specific 
models of f{R,G) gravity were constructed to account 
for the late-time cosmic acceleration p^, and the 
respective constraints of the parameters of the models 
were also analyzed in • So, in order to proceed along 
with the interests of these models, in the present work 
we shall consider two forms of f{R,G) introduced in 
[l0| . In Ref.[l3] the authors have studied the stability 
of de Sitter and power-law solutions in /(i?, G) gravity 
and have shown that gravitational action plays a very 
important role in the stability of the solutions both of 
which depending on the form of the f{R, G) theory 
and the parameters of the model. In this context, we 
further consider the constraints imposed by the energy 
conditions and verify whether the parameter range of 
the proposed models considered in are consistent 
with the energy conditions. More specifically, we define 
generalized energy conditions for f{R, G) modified 
theories of gravity, and consider their realization for 
flat Friedmann cosmological models. In particular, we 
analyze whether the weak energy condition is satisfied 
by particular choices of f{R, G) which were advocated 
inRefs. [3. 

In Sectionini we introduce the gravitational field equa- 
tions for the f{R,G) gravity. In Section ITTll we obtain 
the inequalities corresponding to the energy conditions. 
In Section HVl we consider two specific forms of f{R, G), 
and analyze the constraints resulting from the energy 
conditions. The paper ends with a brief conclusions in 
Section |Vl We use the units c = Gm = 1- 
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II. FIELD EQUATIONS OF f{R, G) MODIFIED 
GRAVITY 

Let us start by writing the most general action for 
modified Gauss-Bonnet gravity 

^^hj d^^V^9f{R, G) + SM{g^\ V') , (2) 

where Suig^^^i^) is the matter action, and f{R,G) is a 
function of the Ricci scalar and Gauss-Bonnet invariant 



defined by 

G = R^~AR^,R^'' + R^,o.pR'""'^. (3) 

Variation of the action ([2]) with respect to the metric 
provides the following gravitational field equation 

Rt.. - y^.uR - n T^^^*) + E^,, (4) 
where S^i, is defined by 



S^. = V^V./fl - g^,afR + 2i?V^V,/G - 2g^,RafG - AR^^x^uIg - AR^^^x^^fo + ^Rt^.^Ia 
+Ag^,R^^V^VpfG + 4i?^a^i.V"V^/G - \ g^,V + (1 - h) (i?^. - ig^.i?). 

I 



(5) 



Note that 



_ df{R,G) 

JR. = ^ 

Or 



_ df{R,G) 
fa = (6) 



V = fiiR+fGG—f{R, G), and Tj™^'-* is the stress energy 
tensor describing the ordinary matter. 
Now, we consider the flat FRW metric 



(7) 



where a{t) is the scale factor. In the FRW background 
with a perfect fluid equation of state for ordinary matter, 
the fleld equations for /(i?, G) gravity are given by 

IrH = -^{p + p) + \{HfR- jR + AH^'fG-mHfG 



-AH^fa), 



(8) 



/rH^ = -p+-{fRR-f-6HfR + GfG-2AH^fG). 



R = 6 i^2H^ + Hj , 
G = 24i/2 (h^ + H 



where p and p are the energy density and pressure of 
ordinary matter, respectively, and the overdot denotes a 
derivative with respect to the time coordinate, t. More- 
over, we have 



(9) 
(10) 



and the gravitational fleld equations may be rewritten in 
the following form 



p,g = -H\ p,s = --(2H + 3hA , (11) 

K K \ / 



where PeS and pos are the effective energy density and 
pressure, respectively, defined by 
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Poff = 


h. 
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Pes = 


Jr 
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Combining the above equations, we obtain the following useful relationship 



(12) 
(13) 



1 

Pes + PeS = -r- 
JR 



P + P+-( -HJr - AH'fG + fR + SHHfG + AH'fc 



which will be used throughout the paper. 
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III. ENERGY CONDITIONS 

In general, the energy conditions arise when one con- 
siders the Raychaudhuri equation given by 

df) 1 

:r = - o - '^M-^''' + ^/^-w''' - Rf^.k^k'' , (14) 

ar z 



where 6 , ct^'^ and are the expansion, shear and ro- 
tation, respectively, associated with the congruences de- 
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fined by the null vector field k^^. 

From Raychaudhury equation it is seen that for spatial 
shear tensor with = a^^cr^^'^ > 0, and for any hyper- 
surface orthogonal congruences, which imposes = 0, 
the condition for attractive gravity, namely dO/dr < 
reduces to R^^k^^k" > 0. In general relativity, using 
the Einstein field equations one can rewrite the above 
condition in terms of the stress-energy tensor given by 
T^jjk^k^ > 0. However, in any other theory of gravity 
such as f{R, G), one should know how to replace in 
terms of T^j,, using the corresponding field equations. 

Equation ^ may be written in the following effective 
form 



gravitational field equations and obtain the energy con- 
ditions as follows 



NEC 



PcS + Poff > , 



DEC 



PcS > and p^s ±PcS >0, 



(18) 



WEC ^ peff > and p^s + Pcff > , (19) 



SEC pcff + 3pcff > and p^s + Pcff > , (20) 



(21) 



where the effective stress-energy tensor is given by 



^/ji(mat) 



(15) 



(16) 



The positivity condition R^^k^k'^ > 0, through the 
modified gravitational field equation (jl5p . provides the 
following form for the null energy condition 



T'^^k^'k" > 0. 



(17) 



Moreover, 



it is plausible to impose the condition 
> for ordinary matter, because it implies 
that the energy density of ordinary matter is positive in 
all local frames of references. 

To deduce the energy conditions in the context of 
f{R, G) modified gravity, we use the modified (effective) 



where the notation NEC, WEC, SEC and DEC stand for 
the null, weak, strong and dominant energy conditions, 
respectively. 

In cosmological context, in addition to the Hubble pa- 
rameter H = d/a, it is appropriate to define the deceler- 
ation, jerk, and snap parameters as 



1 d 



1 'a 1 'a' 



respectively; in terms of which, we may consider the fol- 
lowing definitions 



H=-H\l + q) , 
H = H^{j +3q + 2) , 
H = H''{s-2j-5q-3) 



(23) 
(24) 
(25) 



respectively. Using these definitions, the energy condi- 
tions pB |) - (PT|) take on the following respective forms 



NEC : Pae+Pcff ^ P + P+^ {i&fRRRH\j -q-2f + &fRRH^{] +s + Uq + V + ii) + 

96 {-(6g3 + 27 q^ + 2lq + 8qj + 9j - s)/gg+ 

24[4(g2 + 2q+ 1)H^ + 2q^ + 7q + j + 4]/ggg} if*) > , (26) 

WEC : = p + ^ [-/(i?, G) + G/rH^I - q) - SG/rrH^j - q ^ 2) (27) 

-2mSlG - {2AfH\2q^ + iq + j)fGG\ > 0, PeS + PeS > 0, 

SEC : peff + 3peft = P + 3p + ^[fiR, G) ~ GfRH^l - q) + 108 fRRRH%j -q-2f + 

WrrHHj + s + Uq + Aq^ + 11) + 24HSfG + 
288iJ«(-6g3 - 23g2 - l5q - 8qj - 7j + s)fGG 

+ {2A){288)H'\2q^ + 3q + jffGGG]>0 , Pes + Pc«>0, (28) 

DEC : PcS-PcS ^ P-P+ IhfiR^ G) - mfRRRH^j ^q-2f + 6fRRH\~3j + s + 25q + Aq^ + 19) - 

2AH\fG - 96H%-6q'' - 15q^ - 3g - 8qj - 3j + 5)/gg 

-{24)i96)H'^{2q^ + 3q + j)^fGGG]>0, Pcff + Pcff>0, Peff>0. (29) 
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IV. CONSTRAINTS ON f{R, G) GRAVITY 

In this section, we consider the viable HR, G) modified 
theories of gravity which were used in to study the 
stabiUty of cosmological solutions. Stability of power-law 
solutions for f{R,G) models were given by expressions 
([50)1 and (inj), bellow. In the cases where no fluids are 
considered (vacuum) it is seen that stability of the cosmo- 
logical solutions can be achieved by appropriate choices 
of the parameters space p^. However, in this section we 
consider the constraints imposed by the energy condi- 
tions and verify the consistency between the parameters 
ranges of the models considered in [l^ and the energy 
conditions in f{R, G) flat Friedman cosmological models 
obtained here as 

In this regard, we consider two classes of viable f{R, G) 
given by [lOt 

h[R,G) ^ ^lRPG\ (30) 
f2{R,G) ^ kiR + k2R''G"' , (31) 

where fci, fc2, '^i ^^'^ 7 constants and n, /3 

are assumed to be positive. The Ricci scalar and Gauss- 
Bonnet invariant, defined in equations ^ and (jlOp in 



terms of the Hubble and deceleration parameters, can 
also be expressed as 

R = 6ij2(l - q) , (32) 
G = -24H\ , (33) 

respectively. 

Since the inequalities (|26 |) -(f29 |) imposed by the energy 
conditions in f{R,G) gravity are so lengthly, for sim- 
plicity we only consider the WEC in the following anal- 
ysis. Moreover, we take the following observed values 
for the deceleration, jerk and snap parameters Tl], T3] : 
qo = -0.81 ± 0.14, jo = 2.161° j^, and so = -0.22t°:2i. 



A. f^{R,G) = fiRi^G^ 

To impose energy conditions on fi{R,G) gravity, for 
simplicity we consider the vacuum, i.e., p = p = 0. 
Hence, the WEC constraints, i.e., pcff > and pcff+Pcff > 
0, are given respectively by 



n[6H^il^q)f[-2AH^q]^ 



>o. 



(34) 



fi[6H^l-q)f[-2AHS]'' 
+ {6q^ + 24g2 + 2lq + 8qj + 9 j - s 



6(1 -g)3 ' ^' 6(1 -g)2 

7(1 - 7) [4(g2 + 2q + l)i?2 + + 7q + J + 4]7(7 - 1)(7 - 2)" 



(35) 



6q2 



242ij4g3 



> 0. 




FIG. 1: Plots of the weak energy condition for the specific form of f\{R,G) — fiR^G^. The left and right plots correspond 
respectively to p^g > and p^s + PeB > 0. The positivity requirement of the weak energy condition is satisfied in the plots for 
the parameter range considered. 

I 



The exact analytical expressions for the parameter tained because the constraints provided by the inequali- 
ranges of the constants (3 and 7 can not be explicitly ob- 
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ties ([34)) - (|35|) are so complex, so we consider specific val- 
ues for some of the parameters. In the authors have 
shown that the vacuum solutions are stable for /3 > 1 
and 7 < |(5 — 4/3). To verify this result from the energy 
condition's point of view, we plot the WEC as a function 
of /3 and 7, as is depicted in Fig. [1] The positivity re- 
quirement of the weak energy condition is satisfied in the 
plots for the parameter range considered, for the specific 
form of /i {R, G) given by equation (pO)) . 

B. f2{R,G) ^k^R + kiWC"" 

In this case, as in the previous one, we consider the 
vacuum, i.e., p = p = 0. This modified gravity model ac- 



counts for the stability of the cosmological solutions for 
late-time cosmic acceleration, and the stability of solu- 
tions depends on the values of the coupling constants fci 
and k2- In fact, for ki > 2(ii|^)'^/^fc2 the perturbations 
grow exponentially, and the de Sitter solution becomes 
unstable. For other values of fci and k2 than mentioned 
above, the perturbations behave as damped oscillations 
decaying to zero, hence the solution becomes stable. To 
verify whether the parameter range of the stability of 
solutions are consistent with the energy conditions, we 
again consider the weak energy condition. 
For the form of /2(i?, G) considered by equation (|3ip . 
the WEC constraints, i.e., pcff > and Pcff + Pcff > 0, are 
given by 



A:2(6i?2(l-g))"(-24i/V 



n + m - 



n{n - - q - 2) 

{i-qr 



m{m — 1) 



(2g2 + 3g-f j) 



>0, 



(36) 



fc2[6i72(l-5)]"[-24i74g]'" 
(6g^ + 24^2 + 2lq + 8gj + 9j - s) 



( IV o, (j-g-2)^ , , (j-f s + 14g + 4g2 + ll) 
n{n — l)(n — 2)— j^r; h n[n — 1) 



6(1 -g)3 ' "^■■^ ^' 6(1 -g)2 

m(l - m) [A{q^ -t- 2g -|- 1)7?^ -|- 2q^ -|- 7g -h j -f 4]m(m - l)(m - 2)' 



6g2 



242^4^3 



(37) 
>0, 



respectively. 

It is seen that the week energy condition dose not de- 
pend on ki and k2- Therefore, it turns out that WEC do 
not interfere with the stability conditions. As in the pre- 
vious example, considering the complex constraints pro- 
vided by the inequalities (1551) - (| 3 7 1) , finding exact analyti- 
cal expressions for the parameter ranges of the constants 
m and n is not an easy task. Hence, we consider specific 
values for the parameters to find WEC condition as a 
function of the parameters m and n. The plots shown in 
Fig. [2] again prove that the specific form /2(i?, G) given 
by equation (|5T|) is consistent with the WEC. 

V. CONCLUSIONS 

In this paper, we have studied the viability of an inter- 
esting alternative gravitational theory, namely, /(i?, G) 
gravity. We have considered two realistic forms of 
/(i?, G), analyzed in the literature, accounting for the 
late-time cosmic acceleration and the stability of the cos- 
mological solutions [lo| . The general inequalities im- 



posed by the energy conditions are obtained. To be spe- 
cific and for simplicity, we have focused on the weak en- 
ergy condition and used the recent estimated values of 
the Hubble, deceleration, jerk and snap parameters. We 
have shown the consistency of the above-mentioned forms 
of /(i?, G) with the weak energy condition. For simplic- 
ity, we have just examined the vacuum case for which 
p = p — 0, although this is not a physically relevant case 
bearing in mind that the universe contains matter and 
radiation. This simplification, however, does not change 
the general results, because we can always add a posi- 
tive energy density or pressure satisfying the WEC (from 
matter and/or radiation) to any model which satisfies it- 
self the WEC. Therefore, if the vacuum model satisfies 
the WEC, /(i?, G) gravity with matter does the same. 
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